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Bimodality as a signal of Liquid-Gas phase transition in nuclei? 
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We use the HIPSE (Heavy-Ion Phase-Space Exploration) Model to discuss the origin of the bi- 
modality in charge asymmetry observed in nuclear reactions around the Fermi energy. We show that 
it may be related to the important angular momentum (spin) transferred into the quasi-projectile 
before secondary decay. As the spin overcomes the critical value, a sudden opening of decay chan- 
nels is induced and leads to a bimodal distribution for the charge asymmetry. In the model, it is 
not assigned to a liquid-gas phase transition but to specific instabilities in nuclei with high spin. 
Therefore, we propose to use these reactions to study instabilities in rotating nuclear droplets. 

PACS numbers: 24.10.-i,25.70.Mn,64. 

Keywords: Heavy-ion induced reactions, Fermi energy, Phase transition, Bimodality 



In recent years, the possibility to observe phase tran- 
sitions in finite, even small systems, has received an in- 
creasing interest 0, 0, HI, • In this context, several con- 
ceptual questions are addressed, as the extrapolation of 
thermodynamical properties, or more generally Statisti- 
cal Mechanics, from finite to infinite systems. For in- 
stance, a system with phase transition has discontinuity 
in the equation of state in the transition region. Such a 
discontinuity is not present in finite systems but is ex- 
pected to be replaced by anomalies in specific statistical 
quantities. One of the possible signature of liquid-gas 
phase transition in finite systems is the appearance of 
abnormal fluctuations of the kinetic energy |5( in the mi- 
crocanonical ensemble, these being at the origin of the 
so-called negative heat capacity. Equivalently, a bimodal 
behaviour, i.e. two "bumps" in the energy distribution 
is expected in the coexistence region, the system being 
treated canonically. Bimodalities in event distributions is 
even sometimes promoted as one of the definition of phase 
transition in finite systems 0, Indeed, it is related 
to the anomalous curvature of the entropy or any rele- 
vant thermodynamical potential depending on the con- 
straints upon the system (see however 0). Moreover, 
several studies have shown that bimodality is rather ro- 
bust with respect to the introduction of additional con- 
straints 8] on the system or the long-range coulomb force 
0. In this context, nuclei appear as possible candidates 
to observe liquid-gas phase transition in finite quantum 
system. Indeed, a large var iety of experimental studies 

U^lulllllllllllllllllllllElllini^l^ reports an 
accumulation of "evidences" of critical signals. Among 
them, bimodality in the charge asymmetry of fragments 
produced in heavy-ion reactions at Fermi energies have 
been recently reported for the Quasi-Projectiles (QP) iso- 
lated in peripheral reactions 0, . This signal is pre- 
sented as one of the most robust evidence for the liquid- 



gas phase transition in nuclei. However, the extraction 
of critical signals from nuclear reactions is far from being 
simple. For example, one expects that a signal initially 
present at the chemical freeze-out (i.e. when nuclei do not 
exchange particles anymore) will be largely distorted, or 
even completely washed out by the secondary decay. This 
raises the fundamental question of the phase-space ex- 
plored during the reaction just after fragment formation 
and its modification due to secondary emission. Recently, 
the phenomenological HIPSE 1 model |23| has been devel- 
oped to address these aspects. In this model, very specific 
randomness hypotheses are retained to form clusters in 
the first instants of the reaction, while information on 
the phase-space explored before and after secondary de- 
cay can be accessed without ambiguity. This model, as 
well as the recently developed version for nucleon-induced 
reactions (called n-IPSE) [24(, has been shown to remark- 
ably reproduce experimental observations. 

In this work, we use the HIPSE model to address the 
question of the origin of the bimodality signal in nuclear 
reactions. First the experimental protocol used in refs. 
[2l],|22 is recalled. It is applied to events generated with 
the HIPSE model showing that bimodality is found. Fi- 
nally, we use the possibility to access the phase-space 
before the secondary decay to understand the origin of 
bimodality in the model. To do so, we have generated 
10 6 heavy ions collisions for the Xe+Sn system at 50 
MeV/nucleon. The full impact parameter distribution, 
ranging from the grazing to the head-on collisions has 
been generated. A complete description of the model as 
well as a discussion of the hypotheses used for cluster 
formation can be found in refs. 113, 13- ^ n or der to get 
results directly comparable to those obtained with the 
INDRA 47T array, we have filtered the events and used 
exactly the same experimental protocol (event sorting) 
as described in refs. |2l],|22]. We first use a completeness 
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1 available at http://lpccaen.in2p3.fr/theorie/theory_lacroix.html 



2 



(a) E° T 2 <60 MeV 

- 


(b)60<E,"<120MeV ^ 

; 


(c) 120<E°j<180 MeV ^ ~ 


(d) 180<E°j<240 MeV 






(e) 240<E°j<300 MeV 


(f) E,">300 MeV 






0.2 0.4 0.6 0.8 1 


0.2 0.4 0.6 0.8 



FIG. 1: Correlation between the largest fragment Z\ and the 
charge asymmetry 772 for different bins of (contour levels 
are in log. scale). 



criterion. Here, 'filtered' events, corresponding to the 
best detection of the QP, are selected (80% of the pro- 
jectile); this ensures an almost complete detection of the 
QP products. Due to the forward detection acceptance 
of INDRA mainly semi-peripheral reactions between 5 
to 8 fm are then retained. Complete QP events are fi- 
nally sorted by using the transverse energy of the light 
charged particles (Z=l,2) coming from the Quasi- Target 
(QT), noted E%. By doing so, we avoid the obvious au- 
tocorrelations between the sorting observable (QT) and 
the considered system (QP). Note that in the experiment 
as well as in the simulation, the QP selection has been 
made by taking fragments with positive center-of-mass 
velocities [22J; this assumption has been checked with 
HIPSE and is indeed correct; selecting fragments coming 
from the true QP source or with positive center-of mass 
velocities leads to the very same results for this analysis. 

In the study of bimodality pH 12^ , QT transverse en- 
ergy is assumed to be indirectly related to the order pa- 
rameter and is presented as a way of realizing a " canoni- 
cal" event sorting. Although the transverse energy is in- 
timately correlated to the centrality of the reaction, the 
latter assumption is, in our opinion, far from being clear 
because of the associated large mixing of impact parame- 
ters. In the following, we will focus on the correlation be- 
tween the largest and the second largest fragment emitted 
in the forward center-of-mass hemisphere. We then de- 
fine the charge asymmetry between the two largest frag- 
ments r\ z = {Zi - Z 2 )/(Z 1 + Z 2 ) H2- Zx and Z 2 are 
respectively the largest and the second largest fragment 
charges. Thus, r\z is close to 1 for a large asymmetry and 
it will be the case if an evaporation residue persists after 
de-excitation. By variance, if r\z ~ 0, it corresponds to a 
symmetric fragmentation. Fig. ^displays the correlation 
between Z\ and r\z for different QT transverse energy in- 
tervals. Fig. Ola) shows a single component located at 



r\z ~ 1 and Z\ m Z pro j. This case corresponds to the 
evaporation residue (ER) of the projectile. In Fig. IHb) 
to Fig. Up), we observe a different component, located 
this time at r\z ~ and Z\ ~ 15; the corresponding 
mean fragment multiplicity is here greater than 2, corre- 
sponding to the multifragmentation regime (MF). In Fig. 
IHc). the correlation clearly exhibits both components, 
the asymmetric case (ER) and the symmetric one (MF). 
This coexistence has been assigned to a bimodality sig- 
nal in the fragmenting nuclear systems |22T |. Indeed, by 
projecting the two-dimensional distribution either on the 
x-axis or y-axis, two bumps are observed respectively in 
the distribution of r\z and Z\ (not shown here) for the se- 
lected intermediate transverse energy. These results are 
similar to those obtained in the experimental case [12, 
where a bimodality in r\z has been reported for the QP 
events. 

Let us now specify the properties of the two event 
classes observed in Fig. ^b) and^c). The charge parti- 
tion may mix contributions from quasi-projectile and/or 
mid-rapidity emissions. Indeed, it results from a possi- 
ble simultaneous emission leading to a fragmented freeze- 
out configuration followed by a sequential decay of each 
primary fragments. In the HIPSE model both effects 
are taken into account and the initial partition before 
secondary decay can be easily identified. In particular, 
the number of sources contributing to the bimodality is 
of major interest. By tracing back the origin of clus- 
ters in HIPSE during the decay, we have observed that 
the bimodality signal displayed in Fig. ^ is dominated 
by events where a single source is formed in the for- 
ward center-of-mass hemisphere, this source correspond- 
ing indeed to the excited QP. To get a deeper insight 
on the origin of bimodality in HIPSE, it is necessary to 
clearly identify the phase-space explored by the QP be- 
fore de-excitation for the two event classes. Here, we 
concentrate on Fig. HJc) and we will refer to ER for 
events with n z > 0.8 and MF for n z < 0.2. Fig. H 
presents respectively from top to bottom the correlation 
between the size, the thermal energy, the transferred an- 
gular momentum and the impact parameter for ER (left) 
and for MF (right). The first remarkable aspect ap- 
pearing in top panels of Fig. |3 is that the source sizes 
are not significantly different between ER (Zqp k 50) 
and MF case (Zqp « 45) and the two distributions 
strongly overlap. By contrast, the thermal energy E t h 
is much higher in the MF (E th /A w AMeV) than in 
the ER case (E t h/A w 1.5MeV). It is worth noticing 
that such a result is at variance with a geometrical sce- 
nario like 'abrasion-ablation' models where such an in- 
crease of thermal energy is accompanied by a strong de- 
crease of the QP size j25j. Indeed, in the HIPSE model, 
while the QP and QT are initially formed using geometri- 
cal arguments, the abrasion picture is partially (or even 
completely) relaxed by allowing nucleon exchange and 
by the stron g re organisation due to Final State Inter- 
action (FSI) [2(j. For the beam energy considered here 
(50A.MeV), an important exchange of particle between 
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FIG. 2: Correlation between the atomic number (top row), 
the thermal energy (middle) and spin (bottom) of the QP 
source for events of Fig. ^c). The left column corresponds 
to the ER case and the right to MF case (contour levels are 
in log. scale). 



the target and projectile should also be accompanied by 
a large transfer of orbital into intrinsic angular momen- 
tum (spin). This is indeed confirmed in bottom part of 
Fig. |21 where the correlation between the QP angular 
momentum (noted J) and the impact parameter is dis- 
played. Again, we observed that while ER corresponds 
in average to J w 30ft, in MF much higher angular mo- 
menta are obtained (J « 80H). Consequently, the two 
classes of events associated to bimodality are issued from 
the de-excitation of a single-source with more or less the 
same mass but rather different initial conditions in term 
of thermal energy and angular momentum. Note that 
here, the de-excitation is performed in both cases using 
the statistical sequential decay model of ref . |53 • There- 
fore, in HIPSE, the appearance of the two contributions 
in Fig. (l]is a direct consequence of the statistical decay 
accounting for the initial properties of the QP. To go fur- 
ther, we have performed statistical simulations for a given 
nucleus (Z = 50) by varying the initial thermal energy 
and spin. Fig. [3] shows the corresponding rjz distribu- 
tions. If the spin is set to zero (top row), we see that the 
higher the thermal energy, the smaller the charge asym- 
metry is, but we never obtain a crossover to small values 
of r\z- The situation is clearly different when the spin 
is changed. We see in Fig. [31 that at low spin (J = Oft 
and 30ft) and high spin (J = 90ft) only a single contribu- 
tion exists respectively corresponding to high and low rjz 
values. In opposite, for spin J = 60 — 70ft, we observe a 
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FIG. 3: r\z distributions for the statistical de-excitation of an 
120 Sn nucleus for different initial thermal energies Eth (from 
left to right) and spins J (from top to bottom). 



sudden transition where both contributions coexist. This 
transition appears to be almost independent on the ther- 
mal energy. Therefore, the bimodality observed in Fig. 
is related to a transition governed by the spin trans- 
ferred in the collision, but is not of the liquid-gas type as 
concluded in p2|. 

The physical origin of this bimodal behaviour in the 
statistical decay can be inferred from the value of the 
spin associated to the transition. Indeed, it corresponds 
approximately to the limit of stability of a nucleus with 
Z 50 against prompt fission [28j- At that point, the 
fission barrier height becomes comparable to the energy 
of the least-bound particle and the nucleus cannot resist 
to the deposited spin. In the statistical model this cor- 
responds to a sudden opening of decay channels leading 
to the low rjz contribution. It is worthwhile to mention 
that the description of such instability through a statis- 
tical model is certainly an approximation. Actually, in a 
complete dynamical description of this instability we do 
expect that the system breaks almost at the same time as 
it is formed. In the HIPSE model, the system is assumed 
to be formed and then explore statistically accessible fi- 
nal configurations; it is certainly a too simplistic picture 
and calls for further theoretical developments. In order to 
check our interpretation, we have replaced the standard 
HIPSE deexcitation stag e by a statistical multifragmen- 
tation model (SMM |2j|, using a freeze-out volume of 3 
times the normal density), i.e. we used the HIPSE phase 
space before decay as an event-by-event input to SMM. 
We present in Fig. 01 the correlation between the largest 
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and the second largest fragment obtained with HIPSE 
(left panel) and the new decay provided by SMM (right 
panel) . They display very similar features indicating that 
both scenarii are compatible with bimodality. Therefore, 
we need additional observables, for instance related to 
the kinematical properties. Indeed, first investigations on 
experimental data indicate that the angular distributions 
of the emitted fragments are changed depending on the 
charge asymmetry selection and this could be attributed 
to different spin deposition into the decaying nucleus [3(| . 
However, angular distributions need to be used very cau- 
tiously because of the difficulty to assess the true origin of 
fragments or particles (mid-rapidity, preequilibrium, sta- 
tistical emissions) ; it could blur the spin effect on the an- 
gular distributions |3l| . Such a study, comparing quanti- 
tatively HIPSE predictions to experimental data for QP 
deexcitation, will be reported in a forthcoming article. 
Nevertheless, following the HIPSE scenario which has 
provided a good reproduction of a large number of ex- 
perimental observations, including angular distributions 
for LCP and fragments |23ll3^ |. we conclude that nuclear 
systems close or beyond their limit of resistance with re- 
spect to spin deposition, can be formed, leading then to 
a bimodal behavior. Therefore, heavy-ion induced reac- 
tions might be a tool to study the emergence of shape 
bifurcation associated to high spin, called the Jacobi se- 
quence |33(. This is a very interesting aspect which has 
not been explored in this context. If this interpretation 
is confirmed experimentally, instead of a liquid-gas phase 
transition, we may have a signature of the so-called Ja- 
cobi transition |34j. which is related to a second-order 
phase transition in the continuous limit and has also its 
equivalent in astrophysical context [3{| . 



FIG. 4: Correlation between the largest (Zi) and the 2nd 
largest fragment (Z2) obtained with the standard HIPSE 
(left) and SMM decay (right). Both simulations are made 
with the same set of (Z,A,E*) given by the preequilibrium 
stage of HIPSE (contour levels are in log. scale) . 
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